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APPEND TX 2a 


CONSISTENCY OF PROBABILITIES AND PROCESS SEQUENCE 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Toronto 


https://archive.org/details/31/761117013508 


CONSISTENCY OF PROBABILITIES AND PROCESS SEQUENCE 


Consider the following processes ordered in the same 
Sequence in the POLSIM model as listed here: "Death", "Birth", 
and "Divorce". We are able to observe only two kinds of 
Statistics, the population at risk, N, at the beginning of 
a time period (e.g. one year) and the frequency of success, 

Sir for each of the processes during the given time period. 

We can construct a number associated with each process, 

r,=5,/N which we term the ratio of the process. Input to each 
process WSeits population at risk (Nn; Nir Or N. as the case 

may be) and output is the altered population due to the process. 
(We cannot of course observe these intermediate populations in 
the real world.) The by-product of the process is the number of 


successes of the process. This is illustrated below: 


N 
iis 
P. Birth, rp teases 
N 
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P Divorce, Yr }-————>5 
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If the initial population is N, then the successes, by definition 


of the ratios of the processes, are given simply by: 


Sees ie (1) 
Ss = NE es 
6. = WE (3) 
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Suppose that in the simulation we use the numbers Pyr P, and 
P.,7 aS shown beside each process as the relevant probability. 
From the "Death" process and the hypothesis that bearing a 
child or becoming divorced does not affect the probability 


of dying, we have, 


Sees (4) 


Z 
i 


and (1-P)N (5) 


By comparison of equations (1) and (4) we have: 


D D (6) 
From the "Birth" process and the hypothesis that divorce does 
not depend on the birth process we have, 


os = ees Cr) 


By comparison of equations (7),(5) and (2) we have; 


P, = r,/(1-r)) (8) 


From the "Divorce" process we have, 


Se PN (9) 


By comparison of equations (9), (6), (5) and (3) we derive, 

Be, = r/(i-rp) (10) 
All of the above equations, (1), (2) p21 , (10) sare consistent 
with each other. Equations (6), (8) and (10) provide numbers 
which, if used as probabilities in our simulation, produce 
unbiased results, given our hypotheses. That is, their 


application as probabilities to the population at risk inputted 


to each process will produce the observed number of successes. 


For example, 
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THE INDETAL YEAR STATE VECros 


THE INITIAL YEAR STATE VECTOR 


Chapter 2 listed the elements of the individual state 
vector, and the data source from which it was derived. This 
Appendix provides additional information. Appendix A.2.1 
lists the complete state vector, along with the definition of 
the individual elements, the code by which each element can 
be interpreted, and the way in which each element was derived 


from the Survey of Consumer Finance. Appendix A.2.2 contains 


the record layout for the actual initial year tape. 
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DESCRIPTION OF THE INITIAL YEAR STATE VECrorR 
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DESCRIPTION OF THE INITIAL YEAR STATE: VECTOR 


lige Previous Year Family Unit Number (LYUNIT) 


(a) Definition 


This is simply an integer defining the family that 
the person belonged to in the previous year. Its 
purpose is to enable the POLSIM model to keep 
track of individuals through time. A person's 
present year family identifier would change, for 
example, if he were to leave home and get married. 
The previous year identifier would then be the 


means by which historical continuity was maintained. 
(b) Code 

The numbers are simply sequential integers. The 

first family on the base tape is assigned the 

1dentifier “one”, the secona “cwo',7 “ete. 


(c) Base Tape Designation 


The number assigned on the base tape is exactly 
the same as the "Family Identifier". 


(see 2 below) 


De Family Unit Identifier AUNIT) 
(a) Definition 


This is a number uniquely defining the family to 


which the individual belongs. 


(b) Code 


The numbers are sequential integers. The first 
family on the base tape is assigned the identifier 


"one", the second "two", and so on. 
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(b) Base Tape Designation 


On the SCF Weighted Work file, a family is uniquely 
defined by the following six variables: the 
primary sampling unit number, he segment number, 
the household number, the multiple household 
number, economic family unit number, and the 

census family unit. All individuals on the SCF 
file who have these numbers in common are assigned 
a common family unit number. Every member of the 
first family is assigned a "one", every member of 


the second a “two", and so on. 
oe Province (PROVIN) 
(a) ‘Definition 


This is the province in which the person resides. 


(b) Code 
1 = Nfld 6 = 0ne- 
Oh eo) Pi atk 7 = Man. 
Ss = Nos. S = Sask. 
4 - N.B. OF ates 
oui Pad) LOs = ube 


(c) Base Tape Designation 


The first digit of the five digit primary sampling 
unit number on the SCF file defines the individual's 


province of residence. 


4. Size of Family (SIZE) 


(a) Definition 


This is the number of persons within a given 


individual's family. 
(b) Code 


An integer number. 
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(Top) 


(a) 


Kc) 


Base Tape Designation 


For every family on the base year tape, this 
number is inferred by counting up the number of 
individuals with a common family unit identifier. 
This number is then assigned to all individuals 


within the given Pamiby:s 


Census Family Relationship (DEPNCY) 


Definition 


This variable defines the person's place in the 
family: whether he is the family head, a spouse 


or a dependent. 


Code 
O= Head 
1 - Spouse 


2 - Dependents 
Base Tape Designation 


The 1971 SCF file defines 6 family relationships: 
head, wife, son or daughter, grandchild, brother 
or sister and brother or sister-in-law, and other 
relative. The latter 5 relations are all with 


respect to the head. 


The base tape maintains the "head" and "wife" 
distinctions, and classifies all children and 
grandchildren as dependents. All records con- 
taining? the: last @ relationships, brother or 
sister and brother or sister-in-law, and other 


relative, are omitted on the POLSIM base tape. 


The number of records involved here is negligible, 


and in no way affects the total population size. 
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Marital Status (MSTAT) 
ee Ae) 
(a) Definition 


This is the person's marital status. 


(b} Code 
ie Sineie 
2 =-Nar Pred 
3 = Other 
4 =- Married but spouse absent 


(c) Base Tape Designation 


The first 3 codes are copied directly from the SCF 
file. The fourth code, "married but spouse absent”, 
is added if the family head's marital status 

indicates that he is married, but upon a search of 


the family, no spouse is found. 
Age (AGE) 


The individual's actual age is coded as an integer. On 
the base tape, age is copied directly from the SCF 
Pole. 


Sex (SEX) 


The individual's sex is coded as: 1 - Male, 2 - Female. 
On the base tape, for adults, the sex code is copied 
directly from the SCF file. Chaldven lees) Tian 4 Wo 
not have a sex designation on the SCF file, “and so it 
must be assigned in the creation of the pase year tape. 


This is done on the basis that the probability of a 
child being a male is .5126. 


Major Source of Income (MAJSIN) 
(a) Definition 
Major source of income is the largest component of 


the individual's income: It 12 inferred divectiy 


from the income variables of the SCF tape. 
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(b) 


Code 

ue Investment Incoe 

ai Farm selfemployment income 

3° Non-farm self-employment income 
4. Wages and salaries 

25 No Income 


Weeks in School (WKSCHL) 


Weeks Employed (WKEMP) 


Weeks Unemployed (WKUNEM) 


Weeks in the Non-Labour Force (WKNLF) 


(a) 


(b) 


Definition 


These variables define the number of weeks during 
the year that the person spent in the four mutually 
exclusive states: school, employment, unemployment 
and non-labour force. The four variables will 


always sum to 52. 

Code 

The actual number of weeks is coded, as an integer. 
Base Year Designation 


The four variables were derived from the SCF 
records of weeks worked, weeks unemployed, labour 


force status, and education level as follows: 


(i) If .a person's labour force status as given an 
the SCF file indicated that the individual 
was in the non-labour force, then the weeks 
in the non-labour force was calculated as the 
weeks in a year minus the sum of weeks employed 


and weeks unemployed. The weeks at school 
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(a) 
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eee) 


(Aas) 


variable was set to zero. If, however, the 
difference (weeks in non-labour force) for 
the above calculation proved to be negative, 
i.e., the weeks employed plus the weeks 
unemployed was greater than 52, then the 
amount that the sum was greater than 52 was 
divided by 2 and this quotient subtracted 
from weeks employed and weeks unemployed. 
(Note must be taken that a slight roundoff 
error will occur.)' Weeks in non-labour force 


in this case, was set to zero. 


If the labour force status, combined with the 
April activity code (see 14 below), indicated 
that the person was in secondary school, then 
the weeks, in school, wasis@—e to 07. pe 
weeks in the non-labour force was computed as 
the remaining weeks in the year (twelve) 
minus the sum of weeks employed and weeks 
unemployed as given in the SCF file. If the 
resulting difference was negative, i.e., if 
the sum of the weeks employed plus the weeks 
unemployed was greater than twelve, then the 
amount by which the sum was greater than 
twelve was subtracted from the weeks in 
school. In other words, it was assumed that 
the individual spent less than 40 weeks in 
school (52-WKEMP-WKUNEM, to be exact). The 
weeks in the non-labour force was set to zero 
while the fields for weeks employed and weeks 
unemployed were left as they were on the SCF 
fise. 


If labour force status and yearly activity 
codes indicated that the person was in post- 
secondary school then the weeks in school was 
set to twenty-eight. The weeks in the non- 
labour force was calculated as the weeks 
remaining in the year (twenty-four), minus 

the sum of weeks employed and weeks unemployed. 


Again, if this value results in a negative 
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number, then the amount by which the sum of 
weeks employed and weeks Ghalptoras was 
greater than twenty-four was subtracted from 
the number of weeks in school; i.e., the 
individual spent more than the allotted 
twenty-eight weeks out of a post-secondary 
school. The weeks in the non-labour force in 
this case were set to zero. The fields for 
weeks employed and weeks unemployed were left 
with their original values as given in the 
SCP file. 


(iv) For a person in the labour force, the number 
of weeks in the non-labour force was computed 
simply as the number of weeks in a year minus 
the sum of weeks employed and weeks unemployed. 
The fields for weeks employed and weeks 
unemployed were left with their original 
values and the value for weeks in school was 


set to zero. 


(v) For an individual whose age was less than 
fourte 1 years the weeks in the non-labour 
force was set to twelve, the weeks in school 
was set to forty and the fields for weeks 
employed and weeks unemployed was set to 


Zero. 
4. Education. (EDUCTN) 
(a) Definition 


This variable defines the highest level of education 


attained by the individual. 


(b) Code 

1 - Grade 9 Pi = Univ 3 

2 - Grade 10 12 - Univ 4 

3 - Grade ll 13 — Univ 5 

4 - Grade 12 14 - Univ 6 

5 - Grade 13 15 = Univ 7 

6G = CAAT. 2 16 - Univ 8 

7 - CAAT 2 17 - Univ 9 

8 - CAAT 3 18 - Univ 10 

9 - Univ ] 19 - Less than Grade 9 
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(c) “Base Year Designation 


The education status of the individual was 
inferred from the education field on the SCF 
file. If the SCF record indicated "no schooling" 
or “less than grade 9", the person was assigned 
code 19, "less than grade 9". If the person's 
labour force status indicated that he was at 
present in school, then his educationsstatus 
became identical to his April Activity SGeatas 
(see below). For all others, the mapping from 


the SCF code to the POLSIM code was as follows: 
SCF POLSIM 


a) High School not complete Grade 10 


b) Vocational HS complete Grade 12 
c) High School complete, 
academic Grade 12 or 13 (Ontario) 
dad) Non-University Some CAAT 
e) Non-University Complete CAAT 3 
£f) University some UNLV of 
g) University Completed - 
Diploma UNIV 3 
h) University Completed - 
Degree UNIV 3 


15. April Activity Status (YRACT) 


(a) This variable indicates a person's activity 
status or "what he is doing" in April of the 
year being simulated. A year, in POLSIM, runs 
from April through March. The reason for this 
variable is to provide continuity from year to 
year, as regards a person's labour force behaviour 


and his progress through school. 
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(c) 


Code 

i Grade 9 13 - Univ 5 

é, Grade 10 14 - Univ 6 

5 Grade 11 15 - Univ 7 

4 Grade 12 16" =" Univ 8 

5 Grade 13 Ie = aay 9 

6 CAAT 1 Lo Ory se) 

vi CAAT 2 19 - Unused 

8 CAAT 3 20> Emir 

2) Univ 21. = Unemp. 
10 Univ 2 22 —- Non=Labour Force 
da Unis 3 23 - Age Less Than 14 
2 Univ 4 


pase Year Designation 


ae 


the person's labour force’ Status on tne sce 


record indicates that he is employed, unemployed, 


Owe 


a member of the non-labour force, then he is 


assigned codes 20, 21 or 22, as the case may be. 


ee 
Abe 
TT 


as 


a) 


he is younger than 14 he is assigned code 23. 
his labour force status indicates that he is 
school, then his place in school is assigned 


ie OLMlOW.Sis 


If the person has not completed secondary 
school, and his age is 17 or less, then he 
is assigned a grade on the assumption that 
14 year olds are in grade 9, 15 year olcs 


ty Grade 10; and SO "On. 


If the person has not completed high school, 
but is 18 or over, then it is assumed he 
is in his Jast year of high sehool (grade 13 


in Ontario; grade 12 elsewhere). 


Those who have partially completed work at 
4 University or ‘community “college are assigned 


CAAT D oF ‘UNIV 2, “as "the case may be. 
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Those who have completed community college 
are assigned Code 8, "CAAT 3". 


e) Those who have completed a diploma or certificate 


program are assigned Code 12, "UNIV 4". 


£) Those who have completed a university degree 


program are assigned Code 14, "UNIV 6". 


16. Weight (WEIGHT) 


(a) 


Debinition 


The base year tape contains the records of approximately 
400,000 individuals, each one representative of 23 

of ‘the Canadian: population.” The weighting ftactor 

is thus simply the integer 50, the number required 

to raise the base year population to the total 


Canadian population. 
Base Year Designation 


The weights given on the SCF records are not equal 
weights. One family may have weight 100, another 
200, and so on. In addition, the SCF weights are 
such as to weigh the sample to 10% of the full 
population, rather than the full population 
itself. 


This kind of weighting system is not consistent 
with: the logic of the POLSIM- model. In ~POLSIM, 
the processes of marriage, birth, divorce, etc., 
require that all persons have equal weights, and 
this in turn requires that the lowest common 
divisor of the given SCF weights be selected as 
the common weight, and that additional records be 
created for any family whose weight exceeded this 
value. For example, a single family of weight 100 
would be equivalent to 20 families of weight 5. 
So 19 additional identical records would have to 


be created on the base year tape. 
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Inspection of the data indicated that the lowest 
"practical" common divisor was 5, and this weight 

was chosen as the common weight’ to be made applicable 
to all individuals on the base tape. (A weight of 

5 to 10% of the total population is equivalent to 

a weight of 50 to the full population.) 


Unfortunately, the true common divisor of all the 
weights was not 5, but one. That is, not all of 
the SCF weights were perfectly divisible by 5. 
This creates problems for the replicating process. 
Suppose that a family's SCF weight was 11. Then, 
strictly speaking, it would be necessary to carry 
2.2 such families, with weight 5, on the base 
tape. Since this is not possible, a random rounding 
procedure was adopted to convert the SCF weights 
to multiples of 5. In the above example, this 
would mean that the original weight would be 
converted to either 10 or 15, with probabilities 
-8 and .2 respectively. Division by 5 would then 
yield either 2 or 3, and this would then‘be the 
number of families added to the equal weighted 


base tape. 


17. Employment Category (TYPE) 
(a) Definitron 


This variable defines the way in which a person 
relates to the labotr foree. eas aso: oF 
catch-all variable that indicates whether a person 
is likely to never become unemployed, is likely to 
be subject to unemployment, is retired and so on. 
The complete list of possible "Type's" is given 


below. 


(by, | €ode 


(i) 14 - This is a "Class A" person (a male who 
is either self-employed, or employed in 
a professional, technical, or managerial 
capacity, and who by definition may 
never become employed) who on retirement 


will receive no private pension. 


+ 


Vlaernias sets 


® bal ih iow 
Kdb64 v a 0 @ Ao, 
a _ ; 
FIMIREd wAXeoes virtt 
el 7 
A — 


Y 
$n 119 \ ia 


28h) 
’ 
wii A. Wo oe) Nm 


ly) 15 = A "Clage 2" person who on retirement 


will receive a private pension. 


(iii) 24 - A "Class B" person (by definition not a 
"Class A" person; someone who is subject 
to unemployment) who on retirement will 


receive no private pension. 


(iv) 25 - A "Class B" person who on retirement 


will receive a private pension. 


(v) 3 - A person who has never been a member of 


the labour force. 


(vi) 4 - A retired person who does not receive a 


private pension. 


(Ved ) 5 - A retired person who is eligible for a 


private pension. 
(vial) 140, 150, 240, 2507, A0yes0. 


These are exactly the same as 14, 15, 24, 25, 4 
and 5 except that the person has not yet had an 
initial income assigned to him, or has just had 
an initial income assigned but has mot yet gone 


through an income transition. 


(c) Base Year Designation 


It is first determined whether a person can be 
assigned a "Class A" or a “Class BB" designatien, 
This is done on the basis of the "Occupation-A”", 
"Occupation-B”,; “and “the Income {1elds on tie Scr 
record. If the individual is a male who is either 
self-employed (as inferred from his sources of 
income), or employed in a managerial or professional 
capacity, then he is assigned "Class A" status. All 
other people with occupations are assigned "Class 

B" status. If a person's occupation field is 

coded "not applicable" he is designated as temporarily 


"unclassified" (Type 3). 
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(a) 


(b) 


(c) 


It is then determined whether the "Class A" and 
"Class B" persons are in pensionable employment. 
This is decided on the basis of a "participation 
rate in pensionable employment" (See the Activity 
Block main text and appendix for details). "Class 
A" persons in pensionable employment are assigned 
"TYPE 15", while those in non-pensionable employment 
are assigned "TYPE 14", Similarly, "Claseun 


persons are assigned TYPE's 25 and 24. 


Unclassified persons are presumed to be out of the 
labour force. If they are receiving pension 

income it is assumed they are retired and "TYPE 5" 
is assigned. If they are not receiving pension 
income, but are 65 or older, they are again assumed 
to be retired. In this case “TYPE 4" as assizqnea. 


If the person is not retired, "TYPE 3" is assigned. 


Employment Income (EMPINC) 


Definition 


This is the total income that a person receives 


from employment. 


Code 


The actual dollar value is coded. 
Base Year Designation 


Employment income is the sum of the values given 
in the "wages and salaries", "military pay and 
allowances", and "net income from self-employment" 


fields of the SCF record. 


Interest and Other Investment Income (INTRST) 


(a) 


This is the total income a person receives in the 
form of interest plus all other investment income 


with the exception of dividends. 
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(b) 


ve) 


Code 

The actual dollar value is coded. 

Base Year Designation 

This variable is the sum of the values given 


in the "interest" and "other investment income" 
fields of the SCF record. 


20. Dividends (DIVDNS) 


21. Retirement Pension, Superannuation, and Annuities (RETIRE) 


22. Other Money Income (OTHER) 


(a) 


(b) 


Ge) 


Pete nikterom 


These are the total incomes a person receives 


in the designated forms. 


Code 


The actual dollar values are coded. 


Base Year Designation 


The relevent fields of the SCF record are copied 


directly. 


Dave Total Income (TOTAL) 


(a) 


(b) 


(c) 


Definition 

This is the total market income a person receives. 
Code 

The actual dollar value is coded. 

Base Year Designation 


The required value is obtained by summing all 


of the incomes listed above. 
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POPULATION TABLES 
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Table A3.1 


Survey of Consumer Finances Population of Canada, April 4) 1972 


Seas | (ee es See | 


{ 


| 


Age Group Male Female Total 
0-4 1,801,010 | 1,801,010 | 
am eri ice 2,491,450 
2,102,890 216,070 2,318,960 
1,046,910 1,026,530 2,073,440 
905,090 912,660 | 1,817,750 
784,440 805,030 | 1,589,470 
655,630 688,190 1,343,820 
632,010 626,580 | 1,258,590 
652, 330 638,660 | 1,290,990 
45-49 598, 830 619,690 | 12169520 
50-54 527, 650 526,480 | 1,054,130 
55-59 449, 880 435,410 | 885,290 
60-64 368, 200 372,290 | 740,490 
65-69 322,050 341,180 | 663, 230 
70-74 213,160 248,120 461,280 
75-79 135, 230 176,200 311,430 
80-84 83,790 94,190 177,980 
85-89 32, 980 40,700 73,680 | 
90-94 | 11, 450 10,620 22070). § 
95+ | 1,960 3 150 5,110 
Total iz 13,516,940 PB digs 21, 298,690 
Source: 1972 Survey of Consumer Finances 


Table A3.2 


1971 Base Year Model Population ot Canada 


Age Group Male Female Total 
0-4 922,900 879,600 | 1,802,500 
§=9 1,119,850 | 1,068, 250 | 2,188,100 
10-14 1,179,850 | 1,138,850 | 2,318,700 
15-19 i O47 7250 | 1,027, 350 | 2,074,600 
20-24 905,550 | 912,650 1,818,200 
| 25-29 784, 400 | 804,700 | 1,589,100 
| 30-34 654,950 | 687,850 | 1,342,800 
| 35-39 | 631,400 | 626,050 | 1,257,450 
| 40-44 | 651,900 | 638,600 | 1,290,500 
| 45-49 | 598,750 | 618,400 | 1,217,150 
50-54 | 527,500 i 526,750 : 1,054,250 
| 55-59 450,550 | 435,250 : 885,800 
60-64 | 368, 750 | 372, 200 | 740,950 
65-69 : 321, 400 : 341,700 | 663,100 
| | 
70-74 | 212, 850 | 247,850 | 460,700 
75-79 135, 300 | 175,950 | 311,250 
80-84 83, 450 i 94,200 177,650 
85-89 33,050 | 40,400 | 73,450 
90-94 11, 200 : 10,700 21,900 | 
| 
954 | 1,950 a 3,100 | Bae | 
Total | 10, 642, 800 | 10,650,400 i 21,293,200 | 
Source: Produced by equal weighting routine from the 1972 Survey of 


Consumer Finances 


md 
hore 
a ae 


5S. (neve® 57) old Geil aaisht 
: ‘ofl ; 


~ 
-~ 


° 4 Whe 
DEY ina? 
EL Te @ 


Teun ML 


ek Aa 
se is 
pone 
digh..)4 
mink A , 
atv 41 


_—-)paa:~<- 
. IMF Gis Ag ied 4 


tite 
| de-08 


i ect! 


— $< — 


=a 


"CRD ay lew oo Matias? fepiTpet 


Se’ SOQ)" ape) 


vl 


Table A3.3 


Source: 


Census Division, 


Statistics Canada 


Age Group Male Female Total 
ose | 929,600 | 886,555 1,816,155 
2, | bl sg 7.05 1,101,580 2,254,005 
10-14 | 1,181,450 | 1,129,285 Ors 303,735 
15-19 | 0740435 | 1,039,915 2,114,350 
20-24 941,775 947,630 1,889,405 
25-29 | 800,710 | 783,410 1,584,120 
30-34 | 660,870 : 644,550 1,305,420 
35599 | 645, 045 | 618,825 1,263,870 
40-44 | 640, 770 | 621,760 1,262,530 
45-49 | 613,415 | 625,630 | 1,239,045 
50-54 | 518,900 | 533,640 1,052,540 
55-59 472,415 | 482,315 954,730 
60-64 | 381,695 | 395, 320 777,015 
65-69 | 296,050 | 323,905 619,955 
70-74 | 205,575 251,800 | 457,375 
75-79 | 139,995 185,515 325,510 
80-84 85,680 118,490 204,170 
85-89 40,625 59,390 100,015 
90-94 11, 355 18,815 30,170 
954 2,585 4,620 7,205 
Total 10,795,370 10,772,950 | 21,568,320 


me tae 
Od, . AT 
moi, oe 
' 
yf Sep: 
mG Cad 
Ds 
iy 
oe 
(Pas om 


| so i | 


Wt) Pez 
Wa he 
LT au 


ce » NCD 


Table A3.4 


Relative Po ati 5 
netative Population Errors of 1972 Surve of Consumer Finances* 


(Percent) 


Age Group 
Male Female Total 


* Relative to the 1971 Census 


Source: Tables A3.1 and Isa 3! 
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Table A3.5 


ion* 


r Model Populat 


71 Base Yea 


Relative Errors of 19 


(Percent) 


Female 


Male 


Age Group 


Total 
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Tables A3.2 and A3.3 


Source: 
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A listing of all of the data used by the Immigration 
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VALIDATION TABLES 


The accompanying tables compare simulated distributions 
With actual distributions for level. of education, number of 
children, income class, and actual numbers of immigrants. All of 


these distributions are stratified in the same way that the relevant 


data set is stratified. The actual values are presented first, with 


the corresponding simulated value given in brackets underneath. 
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COMPUTER PROGRAM IMMIG 
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APPENDIX C 


THE DEMOGRAPHIC BLOCK 
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APPENDIX C.1 


THE DEMOGRAPHIC BLOCK PARAMETERS 
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THE DEMOGRAPHIC BLOCK PARAMETERS 
Death Process Probabilities 


The underlying probabilities for this process are 


the survival probabilities: 


oS =" 5 (d,s; t) = Cla,s)+y(a;s) expi=-A (a;,s).(t=1970) > (1) 
as antrodquced in section, 4.2.1, where a=0,1, 2.06... noo 

is age, s is sex and t is time indicating the exact 

year in question. The Meda ea parameters C(a,s), 

Vila, S) gand Ala,s), are given, in Table C.l. im this 

table each combination of sex and age occupies two 

lines. In the first line are shown sex, age, the seven 
estimated survival ratios by W. Zayachkowski, the 
constant C, “the constant y and ‘finally the constant a 

as estimated by our model described in section 4.2.1. 

The second line contains the seven computed survival 
ratios by the formula (1) above. Underneath the y 
constant is shown the t value which indicates the 
significance of the parameters. Finally, the correlation 
coefficient is presented. For degrees of freedom being 
equal to 7-1=6 and at confidence levels of 90% and 95% 

we Hane eritical values of 1.943 and 2.447, respectively. 
As can be seen, all our t-values are significantly 


higher than these critical values. 


For the perjod ‘cf S years rrom 1969 to 1972 the 
survival probabilities are computed and listed in Table 
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Emigration Process Probabilities 


The method used to estimate the emigration 
Probabilities was outlined in section 4.2.2. Te wild 
be recalled that the probabilities are estimated from 
statistics for the period April) 1, 1966 tc Aprel sis) loco 
ang consist of ratios of emigrating individuals €o total. 
population. These ratios are stratified by age, sex, 
and marital) status, and are presented in Table C.3., 


which follows. 
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Birth Process Probabilities 
eee ee  POREDT LILES 


These probabilities are provided by L. Stone, and 
are listed in Table C.4. They are given by single 
years of age from 14 to 48 (inclusive), by birth order, 


and by legitimate-illegitimate cases. 
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Divorce Process Parameters 
eee ew SE Eee Cer SS 


These probabilities are again supplied by 
L. Stone and are given by single years of age> (Ls. tO 59) 
and sex. However, for 60 years of age and over we were 
able to obtain some estimates of divorce rates from 
1968 unpublished statistics of the Vital Statistics 
Division of Statistics Canada. These enabled us to 
enlarge our age dimension up to age 70 and over. 


Table C.5. contains these probabilities. 
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Marriage Process I Parameters 


The parameters of this process consist of probabilities 
that an individual of certain characteristics will get 


married within a period of a year. 


anese probabilities are estimated by five year age 
brackets, sex, region, and marital status (ie, 7 Ssangle 
or "other"). The original data was obtained from Vital 
Statistics Division of Statistics Canada for the year 


1971. Table C.6 which follows contains these probabilities. 
Marriage Process II Parameters 


The parameters’ of this process are used to match 
"Marriageable" individuals into couples. These fall 
inGoetwo séts.. The first gives the probability that an 
individual will choose a mate in age bracket a = 15-19, 
20-247, 6. -a¢rm CONditional on his/her sex and age bracket 
a eew bo 19 w 20-2406 5. e. 5 - The second gives the probability 
that an individual will choose a mate of education e = 

under Grade 9, Grade.9 - Grade 13, College or University, 
conditional on his/her sex and education. These two 
setsvare Jisted in Tables C.7.a and C.7.b., respectively. 
The first set was estimated from Vital Statistics 
unpublished data for 1971 while the second set was 
estimated from the distribution of already existing 
married couples in the Consumer Finance Survey Tape for 
the-vear 1967. It is*fair to say that since the education 
probability distributions were computed by using the 
existing stock of married couples, rather than flows, 
i,@-,enewemarriages within a period of a year, they 
constitute an educated guess, rather than an estimate 


of the desired probabilities. 
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Family Independence Probabilities 


The estimation procedure used to seen s the 
probabilities that a dependent individual will become 
independent and leave home is described in section 4.2.7, 
The estimation is based on inferences made from a Simple 
model in which population is Kept -consteii, st 4 
which mortality, immigration, and emigration are ignored. 
Table C.8 contains the results of the analysis presented 
in section 4.2.7. To increase accuracy, the counts of 
individuals were grouped into two year age brackets as 
Shown in the table. The consequence of this is that the 
probabilities computed by means of formula (EZ) Or 
section 4.2.7., are the probabilities that individuals 
will become independent within a period of two years. 
These are called A, probabilities. If A, are the 
corresponding probabilities that individuals will become 
independent within a period of one year it is easy to 


verify that: 


(1-A,)? = 1-A, or, A, = 1 - V1-A, 


ineseswvesults are given in the last column of eablelc.s. 
Finally, aS waS mentioned in section 4.2.7., our model 

failed to provide reasonable estimates for ages over 26. 
For this reason a probability of 5% is assumed for both 


sexes in the 26 and over age group. 
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TABLE C.8. 


1,908 
4,536 
5,671 
6,033 


ar 67 


1,456 
4,742 
5,028 
3,470 


2,090 


Source of Data: 


Independent Dependent 


36,826 
29,886 
15,663 
pee | 


3,579 


32,246 
20,384 
10,125 

4,571 


2,555 


Married 


340 
2,488 
gf VANS) 

17,804 


20 fa s2 


Independent Dependent Married 


Db a5 
8,345 
19,342 
25,648 


25,991 


ESTIMATES OF FAMILY INDEPENDENCE PROBABILITIES 


Consumer Finance 1969 Survey Individuals' Tape. 
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Interprovincial Migration Process Parameters 


As explained in section 4.2.8 the probabilities 
describing the interprovincial migration process 


consist of: 


(a) Brobabilities q)(a,;ijr) that an individual on 
age a, income i, and who resides in region r will 
stay in his present region of residence. These 


probabilities are given in table C.8va.; 


(ii) Probabilities Ga (ree) that an individual who 
resides in region r will move to region xr" conditional 
that he/she moves out of his/her present region of 
residence.| This conditional transition matrix 1s 


WaLSeacl suis wR CaBalos 7 


(iii) Probabilities G3 (Pp. 9P) that an individual residing 
in province p will move to the Atlantic province 
bay Conditional on the fact that he/she moves to the 
Maritimes. These conditional transition probabilities 


are Lasted| 1n table Cis vc; e 


(iv) Probabilities d4 (PyrP) which are the same as (i111) 
above for the Prairies provinces' transition. These 


probabilities are listed in table C.8.d. 


Clearly, the above sets of probabilities imply the 
existence of transition matrices from province to province 
which are conditional on income and age of the individuals. 


These transition matrices are listed as tables C.9.l, 
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Validation of Demographic Block Parameters = 


Section 4.3.4 outlined the general way in which 
the parameters of the Demographic Block could be validated 
and adjusted. It will be recalled that the method 
entarled, first of all “the abeearhene. fe of a confidence 
interval for the "successes" of a given process. This 
confidence interval depended on the observed population 
at risk and the probability of success for the event in 
question. If the observed number of successes of a 
given process was then within the relevant confidence 
interval, the probability of success could be deemed 
"valid". If not, it would have to be adjusted by a 
correction factor, which was a function of the observed 
number of successes, the population at risk, and the 


Possibility of success, 


The present appendix presents the details of this 
procedure for the processes of the Demographic Block. 
For any given process, tables of populations at risk, 
observed successes, probabilities of success, and the 
implied confidence intervals are presented. The correction 
factors that are derived from these parameters are then 


given, along with tables of the new adjusted probabilities. 


In some cases, such as the survival process, the 
above procedure is not feasible. In these instances a 


more aggregate validation is presented. 
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Survival Process Probabilities 
Bee OAD L LACS) 


The Survival Probabilities for almost all ages are 
very close to unity. The method of validation described 
in section 4.3.4. for these processes failed to produce 
any fruitful conclusion for two reasons. First, the 
initial population errors are relatively quite high; and 
second, the probabilities either of death or of survival 
are very close to zero or to unity, respectively. Both 
reasons made the probabilities very sensitive to any 
adjustment, and it was found that either the death 
probabilities required a high order of correction or that 


the complementary survival probabilities exceeded unity. 


However, a validation was made at the most aggregated 


level and it produced the following results: 


For the year 1968 we expect 80,239 deaths for males 
and 55,676 deaths for females out of populations of 
9,946,150 males and 9,951,850 females. These totals 
result from an analytical application of the estimated 
probabilities for the year 1968 to the 1968 Consumer 
Finance Tape population. The death ratios for this year 
are .0081 for males and .0056 for females. In the same 
year the Vital Statistics figures are .0086 for males and 
-0062 for females. This indicates that our probabilities 
are optimistic and on the aggregate we will loose 5 to 6 
deaths for every 10,000 of population. This error is 
within our tolerance set by the simulation error margins. 
Therefore, despite our failure to follow a more comprehensive 
procedure of validation, the aggregate picture is quite 


satisfactory. 
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ee 2s Emigration Process Probabilities * 


With these prenaniiieres we had the following 
conceptual problem. From statistics of emigration flow 
for individuals we derived probabilities on Snqration 
of families based on the characteristics of the family 
head. Comprehensive validation of such estimates is not 
possible due to the absence of statistics of emigration 
flows of families. However, we analytically estimated 
that for the year 1968 we should expect a total of 
67,883 emigrants, while the actual aggregate emigration 
was recorded to be 67,270. The error is of the order 
of 0.9% (overestimate) which could be judged as satisfactory. 
The estimated joint distribution of emigrants over ages, 
regions, and family sizes might be less accurate, however. 
At the moment, due to the lack of available statistics, it 
is impossible to assess the extent of the errors at the 


more disaggregated levels. 
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Validation of Birth Parameters 


In table C.10. are listed the 1968 female population 
eligible for birth by age and marital status, together 
with the observed number of births for the same year. 
The next table, C.1l., contains the average birth 
probabilities in each 5-year age bracket for each 
marital status. The suggested procedure of validation 
in section 4.3.4. indicated that! alll the birth probans licics 
should be adjusted by the correction factors listed in 
table C.1l. The adjusted probabilities are listed in 
table C.12. From inspection of the table of correction 
factors it seems that L. Stone overestimates the probabilities 
for the legitimate (married females) cases while underestimating 


those for the illegitimate (non-married females) ones. 
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Validation of Marriage Probabilities 


The tables that follow describe the validation of 


the probabilities that an individual will get married 


Within a given year. Table C.13. lists the 126¢ population 


of individuals eligible for marriage by 5-year age 
groups, by sex, and by marital status. In table C.14. 
the observed marriages for 1968 are listed, stratified 
on the same characteristics. Both tables C.13. and 
C.14. are derived from Vital Statistics data. In table 
C.15. the average probabilities are listed while in 
table C.16. the correction factors are presented. In 
table C.17. the confidence intervals are shown as 
Suttaned in section 4.3.4. Finally,eime tabic Coleg 


the adjusted probabilities for marriage are listed. 
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Validation of Divorce Probabilities 

As in the marriage process we list the observed 
population at risk in 1968 in table C.19. by sex and 
beyear age groups. In table C.20. are listed the 
observed number of divorces during 1968. In table 
C.21. the average probabilities are listed, while 
tale ©.22. lists* the corrections factors. | Table: co: 
then presents the confidence intervals, and Fine Liv 


table C.24. lists the adjusted divorce probabilities. 
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ae ae Validation of Interprovincial Migration Probabilities 


The Interprovincial Migeetice Probabilities were 
validated in a different way than were the probabilities 
describing the other Demographic Block processes. The 
reason was the lack of a second source of information. 
im Caple.C.25. there is shown the distribution of 
population by age and province in the year 1968 as 
recorded in the Consumer Finance 1968 tape (April 1968 
population). By using the interprovincial migration 
probabilities and ignoring death, emigration, births, 
and immigration, the new distributions of population 
were computed for the years 1969, 1970, 1971 and 1972 


as=shown in tables C.26, C.27, €C.28 and GC. 297 erespecrively. 


Again from the Consumer Finance 1972 tape (April 
1972 population) we obtain the distribution of population 
by age and province and this is shown in table C.30. 
By comparison of tables C.29. and C.30. we can see that 
the differences between these distributions are not 
significant. This is especially true when account is 
taken that many processes influencing these distributions 
are ignored in the preparation of Table C.29. The 
Interprovincial migration process can therefore be 


deemed satisfactory. 
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APPENDIX C.3 


INITIAL POPULATION ERRORS 
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7 


21h eA PARE 


INITIAL POPULATION ERRORS 


In this Appendix we report initial population 
errors, for the 1968 year, as they pertain to the following 
demographic block processes: the death process, the birth 


process, the marriage process, and the divorce process. 


Table"C.31. contains initial population errors ror 
the whole national population by age, sex and region. Since 
the whole national population is the population at risk for 
death, these errors signify the tolerances on tie number of 
simulated deaths. For example, the 1967 Consumer Finance 
Tape contains 10.8% too few male individuals from Quebec in 
the age bracket 10-14 years of age (relative to Census 
Statistics). Therefore we are bound to have an error in 


this group for deaths of at least 10.8%. 


Table C.32. contains the initial population errors 
for the female population from 15 to 49 years of age. This 
population is the population at risk for giving birth. The 
most significant error in this group is the married population 
of age 25-29 years (-10.2%). The other errors in this group 
which are large, 24.2% and 10.1% for the married population 
of ages 15-19 and 40-44 years of age, respectively, are not 
too significant since the birth probabilities in these 


groups are not as high as for the 25-29 age group. 


soni Prion vet io e 
 gabateine, Sexe! 7) bina 


46 ataris Moisslie=s® fait) wibiins Ae vita 
pow (ed tira? ‘hn SS0 nha «Ve eeritonta banouten sands it J re 


7 S64 ASI, 0) adloni vqow ad2 er Ath among Foti aed 2 
Te sGrlgon' ® Lone ulgsé tec “nee Cbbheds ont ane wept - pi 
urinilot, sagtaey Chel $s (aimee see . assed a > 
Hi ideo’ “att Stavlhylbi-oj/ an wet = on.1 nud eFreG eqat - 
hupiss of e¢isnie:) evel Te mhle bi=0) Se oe ” 
to - Gi COur WS (Neo Sie we wy gaat . tool emeinae ‘ & 
‘UOC soni. t4) ty Cab Ret at qwesy widz 7 
iiws dokfnisyag Nels ind oft aekeehe eke aaa 
Stith: .O06 Io atagy. Ve Gy bl] tea i ala wlnees ert 303 
aiY .ilasid vis vid. te? Ande 3 » WO) elon one o2 iokthtis qo 
old ntarnjicm beks Sum sic wor” ehhi wh Jane srmeitingts nou 
nein. Viel pet B¢ THe Byes - {0G ela) ato'y i¢-"e one Ya 
; PORT ee Toy Sian ali! ee) © oe 0648S epee: ene) Ole 
4% Su 5 esue bisehceka) yap Go Hiney bE-Od bra O12) tone ip 
‘nents ih amv. aideray asd ans eonke simpli Dingla aos 
qo Ty ope, di~eS aly 2A as ag ait en son sin equows 
- . 
; : 


Table C.33. contains initial population errors for 
the married population, which is the population at risk for 
divorce. As can be seen, the Consumer Finance file underestimates 
the number of married individuals for almost all ages and 
for both sexes, except for the age bracket 60-64. A few of 


these errors are quite significant. 


Pinably, stable C.34. contains the initial pepulatiton 
errors for the non-married, which constitutes the population 
at risk for marriage. The errors for single people are 
almost» all negative, i.e., the Consumer Finance file under- 
represents the population of single people. The errors are 
almost all positive for "other" (divorced, widowed, separated, 
etc.), i.e., the Consumer Finance file over-represents the 


"others". 
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